Abstract 

We present a preliminary analysis of time series photometry on the 
pulsating subdwarf B star PG1336-018 acquired as part of two Whole 
Earth Telescope runs. PG1336-018 is a multi-mode pulsator with at 
least 15 periods centered around 175 seconds. It is also an HW Vir-type 
eclipsing binary with a period of 2.4 hours. 
, Though over 20 pulsating subdwarf B stars have been discovered, 

' only a few have resolved temporal spectra and still fewer have identified 

i pulsation modes. We analyze features unique to the eclipsing binary 

Q ■ system of PG1336-018 as a means of identifying pulsation modes. 

1 Introduction 

Pulsating sdB (known as sdBV or EC14026) stars offer us the potential to 
^ ■ probe the interior of Blue Horizontal Branch (or sdB for field) stars using 

. asteroseismology. PG1336-018 (hereafter PG 1336) was discovered to be a 

^ I pulsating sdB star by Kilkenny et al. (1998), with the added bonus that 

it is also an eclipsing binary. The companion is a late-type ~M5 dwarf, so 
the stellar radii are comparable, yet the companion contributes little to the 
integrated flux. During primary eclipse the companion covers about half of 
^ _ the pulsator. 

P ; 2 Previous studies 

i PG 1336 was the target of a Whole Earth Telescope (WET) campaign in 

April of 1999. It was observed to pulsate in at least 15 modes (Reed et al., 
2000) with 2 modes split by twice the orbital frequency. We wished to use 
the eclipse to identify some of the remaining modes. 
I Pulsations manifest themselves as alternating hot and cool regions on 

the star's surface. As the visible surface area changes (during eclipse), so 
should the observed pulsation amplitudes. An analysis in which the pulsa- 
tion axis was aligned with the rotation axis was presented in Reed, Kawaler, 
& Kleinman (2001). They determined that modes with £ = 1, m = 0; £ = 2, 
m = ±1; and all £ = 3 modes would only be visible during eclipse while 
i = 2, m = ±2 modes would disappear during eclipse. Other modes would 
remain unaffected. This means eclipses could be used to identify pulsation 
modes. 

Unfortunately, Xcov 17 did not have the usual WET coverage and was 
not optimized to study the eclipses. However it did prove the method to be 
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viable and, in April, 2001, the WET observed PG 1336 again. 

3 Xcov 21 

During Xcov 21, ~440 hours of data were collected from April 17-30, 2001 
(a duty cycle of '^62%). Over 27 separate modes were identified in the out- 
of-eclipse (OoE) temporal spectrum and this time the data were optimized 
to examine the pulsations through eclipse. Though still a non-trivial task, 
wc were able to detect 15 peaks. As summarized in §2, we would expect to 
observe several modes in both the OoE and primary eclipse data, as well as 
new, rotationally split modes in the primary eclipse data. This was not the 
case. Of the 15 modes detected during primary eclipse, only 2 were observed 
in the OoE data (the 2 with the largest amplitudes). 

Our first assumption was that the primary eclipse data was too contam- 
inated with aliases to be useful. To test this, we reduced and analyzed data 
from the secondary eclipse in the same fashion as from the primary eclipse. 
We readily detected 12 peaks on the temporal spectrum, 11 of which co- 
incided exactly (to within the errors) to those detected in the OoE data. 
This gave us confidence that it was not a reduction problem. It more likely 
indicates that aligning the pulsation axis with the rotation axis is a poor 
assumption. 

4 A tipped pulsation axis 

As first suggested by Margarida Cunha (private communication) , if the tidal 
force dominates the Coriolis force, the pulsation axis may align along the 
axis of the binary. This is similar to the roAp stars where the magnetic field 
dictates the alignment of the pulsation axis. To determine the effect of such 
an alignment, we created simulated light curves where the pulsation axis 
is tipped by 90 degrees (the inclination of the actual binary is 81 degrees) 
and precesses to always point toward the companion (eclipses were ignored) . 
Our simulated light curves and Fourier transforms are provided in Figure 1. 
Pulsation modes are indicated on the left of each plot (as in Urn). For 
the light curves, phases (where a shift of 0.5 = 180") are dashed lines, 
with scales on the right of the plot and dotted lines indicate orbital phases 
(PE=primary eclipse, SE=secondary eclipse, T and A=quadrature when 
the star is moving toward and away from the observer, respectively). The 
dot-dashed line indicates the input frequency (corresponding to a pulsation 
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Figure 1: Left: Simulated light curves and FTs with the pulsation axis 
tipped 90 degrees to the rotation axis. 

period of 175 seconds) of the FT, the dotted line is 1 • forb away from the 
input signal and the dashed line is 2 • forb away from the input signal. As 
indicated in the figure, every non-radial mode undergoes at least two 180" 
shifts in phase during each orbital cycle. 

Each pulsation frequency produces a complex light curve which creates 
a many-peaked temporal spectrum. But this complexity may be the key 
to understanding the pulsations. PG 1336 has many peaks in the temporal 
spectrum; by tipping the pulsation axis, we can create many peaks from 
fewer actual pulsations. If the lightcurve is binned by orbital phase for a 
particular mode, then only the real peak will be present, and predictable 
phase shifts for similar valued modes will be observed. 

Of the 4 highest amplitude modes, two are split by 2 • forb, a signature 
of £=1 modes. We divided up the light curve and binned by orbital phase 
appropriate for an 1=1, m = ±1 mode. We then analizcd each temporal 
spectrum independently. In both (phase) CctSGS, Si peak midway between 
appeared where none previously existed in the combined spectrum. Addi- 
tionally, the phase between the regional spectra changed, as expected for an 
£=1, m. = ±1 mode. This provides strong indications that indeed PG 1336 
has a tipped pulsation axis. The next step will be to bin the light curve by 
orbital phase appropriate for other modes and perform similar tests. Such 
tests will be performed in the near future. 
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Figure 2: Frequencies detected in the temporal spectra of PG 1336 binned in 
amplitude (for clarity). Solid lines are modes present in the OoE temporal 
spectrum while dashed lines indicate peaks only observed in the labeled 
temporal spectra. 

4.1 Eclipses revisited 

Simulations like those in Reed, Kawalcr, & Kleinman (2001) were created to 
model light curves during primary eclipse for the re-aligned pulsation axis. 
With the exception of the £=2, m = ±2 modes, all modes show no effect 
from the eclipse in phase and only mild changes in amplitude. Thus, if the 
pulsation axis is tipped 90 degrees to the rotation axis, amplitude during 
the primary eclipse is no longer an indicator of pulsation modes. 

5 Conclusions and future work 

Figure 3 shows the pulsations detected from various reductions of the data 
for PG 1336 from the Xcov 21 WET run. More than 20 individual pulsa- 
tions are detected in the combined data set, most of which are detected in 
other, regional data sets as well. The primary eclipse temporal spectrum is 
the most incongruent, with only 2 modes detected in other temporal spec- 
tra. Many of the remaining modes are orbital separations away from modes 
detected in the OoE spectrum. 

The temporal spectra at eclipse and quadrature provide good indications 
that the pulsation axis of PG 1336 is not aligned with the rotation axis, 
but rather (at least somewhat) with the tidal force from the companion. 
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This likely has far broader applications than just PG 1336. Since a sizable 
fraction of sdB stars are likely in binaries (Green, Liebert, Sz Saffer, 2001), it 
is probable (and almost a certainty for KPD 1930+2752) that other pulsators 
also have pulsation axes which (somewhat) align with the tidal force of their 
companions. 

I would also like to call on theorists to help in understanding how pul- 
sations respond to tidal forces. It is already believed that p and g modes 
will "feel" the tidal forces differently, but it would be helpful for observers 
to have some form of parameterization so we can predict what the pulsation 
axis should be. There is also some question (at least in my mind) as to what 
it means to have m modes which are not due to rotation. 
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